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This paper describes the development of a 
hybrid computer simulation of a TF34-GE-100 turbo- 
fan engine with post-stall capability. The simu- 
lation operates in real-time and will be used to 
test and evaluate stall recovery control modes for 
this engine. The simulation calculations are per- 
formed by an analog computer with a peripheral 
multivariable function generation unit used for 
computing bivariate functions. Tabular listinqs 
of simulation variables are obtained by inter- 
facing to a digital computer and using a custom 
software package for data collection and display. 
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Nomenclature 


area 

nozzle flow coefficient 
nozzle flow parameter 
specific heat, constant volume 
force 

gravitational conversion constant 

heating value fuel 

enthalpy 

enthalpy drop 

energy 

rotor inertia 

mechanical equivalent of heat 

constant 

length 

rotor speed 

pressure 

pressure drop 

pressure ratio 

gas constant of air 

flow resistance 

temperature 

downstream temperature 
temperature ratio 
turbine flow parameter 
turbine enthalpy parameter 
volume 
stored mass 

flow 

total pressure relative to S-L 
specific heat ratio 
efficiency 

total temperature relative to STD 


Subscripts 


BLD bleed 

C compressor 

CB combustor 

FC fan or compressor 

PH fan hub 

FT fan tip 

I inlet 

N nozzle 


TB turbine 
f fuel 

i initial condition 
in into volume 
out out of volume 


Introduction 

The desire for increased performance from 
aircraft gas turbine engines has resulted in 
engine operation close to stall. Because of this, 
the possibility of engine stall has increased. In 
most instances of stall, there is a momentary 
change in the engine airflow pattern and a loss in 
thrust. Usually, the engine returns to its normal 
operating condition without the need for cor- 
rective control action. This is a recoverable 
| stall. However, in some cases, the engine fails 
to recover and settles at a stable, low-thrust, 
operating point with low rotational speeds and 
high turbine temperatures. To restore power and 
prevent serious engine damage, the engine must be 
shut down and restarted. This is termed nonrecov- 
erable or "hung" stall. 

The reasons why some stalls are recoverable 
and others degenerate into nonrecoverable stalls 
are not fully known or understood. The NASA Lewis 
Research Center is conducting a research program 
to gain a better understanding of stall occurrence 
and to design and evaluate control strategies for 
stall recovery. 

A necessary tool in the understanding of stall 
phenomena is a mathematical model of the engine 
dynamics. Such a model must accurately describe 
both the normal out-of-stall and the in-stall 
behavior of the engine. Previous work has been 
done to develop a post-stall model /simulation 
of the General Electric TF34-GE-100 engine 
compression system. That simulation is capable of 
producing recoverable and nonrecoverable stalls, 
depending on the value of certain engine para- 
meters. However, that simulation is limited to 
operation at fixed compressor and fan speeds since 
the engine turbine models and rotor dynamics are 
not included. 

To evaluate control strategies for stall 
recovery, a simulation that can represent steady- 
state and dynamic engine performance over a full 
range of operating conditions (i.e., out-of-stall 
and in-stall) is required. Furthermore, if the 
simulation is to be used to validate actual 
control system software and/or hardware, then the 
simulation must operate (i.e., solve the dynamic 
equations) in real time. Hybrid computer simu- 
lations have been shown to provide this real-time 
response while retaininq a high dearee of steady- 
state accuracy. 4 


This paper describes the formulation of a 
mathematical model and the development of a full 
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range, real-time, hybrid computer simulation of 
the TF34-GE-100 turbofan engine. 

The mathematical model of the out-of-stall 
engine performance was derived from the engine 
manufacturer's steady-state digital simulation of 
the engine. The in-stall model was based on the 
results of Greitzer's work. 3 - 5 Many of the as- 
sumptions and techniques used to develop the over- 
all hybrid simulation were based on earlier 
developments of real-time hybrid computer simu- 
lations of the Pratt and Whitney TF30-P-3 and 
F100-PW-100 augmented turbofan engines. 6-8 

The engine model was programmed to run on one 
of the NASA Lewis Research Center's two hybrid 
(analog and digital) computer systems. The bulk 
of the simulation calculations are performed on 
the analog computers. A peripheral multivariable 
function generation unit is used to evaluate the 
model's bivariate functions. The digital computer 
was interfaced through analog-to-digi tal con- 
verters (ADC's) and a custom software package was 
developed for data collection and display. 

Engine Description 

The TF34-GE-100 turbofan engine is designed 
and manufactured by the Aircraft Engine Group of 
the General Electric Company in West Lynn, 
Massachusetts. The TF34-GE-100 is a dual rotor, 
front fan, two stream, separated flow, high bypass 
ratio turbofan engine. A single inlet supplies 
airflow to a single stage fan. After leaving the 
fan, this airflow separates into two streams, one 
through the engine core; the other through the 
bypass duct. Since the TF34 has a high bypass 
ratio, the bulk of the airflow passes through the 
bypass duct and when exhausted, provides the major 
portion of engine thrust. 

The core airflow passes through the com- 
pressor which has fourteen stages and variable 
stators. In the combustor, the air is mixed with 
fuel and burned, providing the hot gases to drive 
the turbines. The gas generator turbine drives 
the compressor. It consists of two air-cooled 
stages and is connected to the compressor by a 
hollow shaft. The fan turbine has four stages and 
drives the fan through a shaft concentric to the 
gas generator shaft. 

The bypass duct airflow and core gas flow 
exhaust through separate nozzles. The exhaust 
nozzle areas are fixed and convergent. The core 
stream convergent nozzle is concentric to the 
bypass stream convergent nozzle. A schematic 
representation of the TF34-GE-100 turbofan engine 
is given in Fig. 1. Engine station identi- 
fications are also given. 

Model 

The mathematical model of the TF34-GE-100 
turbofan engine must represent steady-state and 
dynamic engine performance over a full range of 
operating conditions (i.e., out-of-stall and in- 
stall). The model reported herein differs from 
previous models 0- ® in that it provides for in- 
stall operation. This model is based, to the 
extent possible, on data obtained from the engine 
manufacturer's steady-state digital simu- 
lation. y In some cases, other sources of 


TF34 engine data were needed when developing 
models of specific components. Figure 2(A) shows 
the computat ional flow in the model for normal 
engine airflow. In Fig. 2(B), component models 
that must account for reverse flow or rich fuel- 
air mixtures are noted. 

For steady-state accuracy, wide range per- 
formance maps are included for the fan and com- 
pressor. To provide transient capability, fluid 
momentum, mass and energy storage and rotor 
inertias are included. Intercomponent volumes are 
assumed at engine locations where gas dynamics are 
considered important or necessary to avoid 
iterative calculations. Dynamic forms of the 
continuity and energy equations and the state 
equation are solved for stored mass, temperature 
and pressure in each volume. The effect of fluid 
momentum on engine dynamics is also considered in 
this simulation. The momentum equation is solved 
for flow in the fan tip, the fan hub and the 
compressor. 

For flow reversal which may occur during 
stall, the source temperatures for the affected 
energy equations are switched from the normal 
upstream temperatures to the temperatures within 
the particular volume. To avoid discontinuities 
in the corrected speed for the fan and compressor, 
the upstream temperature value is still used as 
the correcting temperature during reverse flow. 
These in-stall models for the inlet, fan, and 
compressor are based on the work reported in Ref. 

1. Fan or compressor stall is detected when, for 
a specified corrected speed, airflow is less than 
the airflow at the stall line. To attain real- 
time response of this simulation and to minimize 
the number of analog components needed, some 
simplification of the basic aerothermodynamic 
relationships in the model was necessary. Fluid 
properties, such as specific heats and gas con- 
stants, were assumed to be constant (i.e., 
independent of temperature and fuel-air ratio) 
over the engine operating range. Additionally, 
the effect of fuel-air ratio on enthalpy in the 
turbine volumes was assumed negligible. 

Descriptions of the individual component 
models follow. General forms of the equations are 
given. 

Inlet 

The inlet is modelled as a pressure drop using 
the general equation form 

*j = aPj/R (1) 

The inlet conditions are set as necessary to 
attain the desired fan inlet conditions. Volume 
dynamics for mass and energy storage are in- 
cluded. The forms for the energy and stored mass 
equations, used in this model, are 

It (WT) - Y * [ Z <*1n * T in> - 1 ("out * T >] (2) 

It (W) " r *irf £ *out (3) 

where ft. denotes flow into the volume, & 

in » out 


2 


denotes flow out of the volume, Ti n is upstream 
temperature, and T is the temperature of the 
volume. The inlet pressure and temperature are 
calculated from the gas state equation 

P*Vs=W*Ra*T (4) 

During flow reversal in stall, the inlet acts 
as an energy sink. As discussed earlier, the 
reference temperatures change and hence, the form 
of the energy equation becomes 

(WT) * y * (i (w in * T) - i (w Qut * T')) (5) 

where T' is the downstream temperature. 

Fan 

Overall performance maps are used to represent 
the out-of-state steady-state performance of the 
fan. There are separate maps for the fan tip and 
hub sections due to the radial pressures gradient 
that exist in the fan. The fan tip and fan hub 
pressure ratios, (P/P )py and (P/P)p H , are 
computed from bivariate functions of corrected 
total fan airflow and corrected fan speed. Cor- 
responding fan discharge temperature ratios, 

(T/T)ft and ( T/T )ph * are calculated as line- 
ar functions of the respective pressure ratios. 

This avoids having to compute fan adiabatic 
efficiency and calculations involving expo- 
nentiation. From the fan tip and fan hub pressure 
and temperature ratios, the fan tip and fan hub 
pressures and temperatures, Ppy, Prh, Trj, 
and Tpn are calculated. Volume dynamics for 
the fan tip and hub volumes are included. The 
equations for the energy and stored mass, 

(WT)pt, Wpj, (WT)ph, and Wpn are similar 
in form to Eqs. (2) and (3) of the inlet. Energy 
contributed by shaft work is accounted for by the 
calculation of the fan discharge temperatures. 

Fan duct and compressor inlet pressure and temper- 
ature conditions are calculated from the gas state 
Eq. (4). Radial crossflow in the engine goose- 
neck, between the fan tip and hub sections is 
assumed to be proportional to the pressure dif- 
ference. Zero crossflow is assumed at the design 
steady-state point. Additionally, flow dynamics 
are considered for the fan tip and fan 

hub flows, ftp-j. and Wp H by use of fluid momentum 

equations. The general form for the momentum 
equation is 


Compressor 

As with the fan, the out-of-stall steady-state 
performance of the compressor is represented bv an 
overall performance map. Inputs to this bivariate 
t map are corrected compressor flow and corrected 
compressor speed. The output of the map is the 
compressor pressure ratio, ( P/P ) r. Compressor 
temperature ratio, (T/T) c , is fit as a function 
of the pressure ratio. Volume dynamics for the 
energy and stored mass, (WT)c and Wc, are 
included. These are similar in form to Eqs. 

(2) and (3) of the Inlet. Compressor flow, , is 

calculated by a momentum equation similar to Eq. 

(6) given in the Fan. 

As for the fan, the in-stall model is based on 
the compressor model contained in Ref. 1. The 
stalled form for the compressor energy equation 
reflects the changes in reference temperatures as 
shown in Eq. (5) for the Inlet. A phi-psi re- 
1 lationship from Ref. 1 calculates the stalled 
pressure and temperature for both reversed and 
forward compressor flow during stall. 

Bleeds 

Compressor discharge bleed is used for turbine 
cooling. The turbine in-flow bleeds are modelled 
simply as constant proportions of compressor 
flow. The general equation form is 


Combustor 

Combustor airflow is calculated directly from 
the pressure drop across the combustor given by 

*CB = V ^CB ^CB^CB * T CB> < 8 > 

An algebraic energy equation replaces the 
dynamic form and is used to calculate combustor 
enthalpy such that 

h CB “ (hw W*CB,out (9) 

where 

(h ”^CB " h CB,in * Vout + *f Mr, * HVF) (10) 


w = (I/O £ (a (P * A) + F) dt + « (6) 

The in-stall fan model is based on the model 
of Ref. 1. The stalled energy equations for the 
fan tip and fan hub volumes provide for changes in 
the reference temperatures as was done for the 
Inlet volume. The form of these energy equations 
is similar to Eq. (5) of the Inlet. 

The stalled pressures and temperatures of the 
fan are calculated from the phi-psi relationship 
given in Ref. 1. This simplifies the fan stall 
model by avoiding bivariate map evaluations during 
stall. Reversed and forward fan flow during stall 
is accounted for in these equations. 


The combustor temperature is computed from a 
bivariate function of enthalpy and fuel-air 
ratio. This function is based on data contained 
in the standard gas tables.^ Since no mixing 
of the flows occurs in the combustor discharge 
volume, a simple first order lag is used in place 
of the dynamic energy equation. The time constant 
is assumed to be equal to the ratio of the stored 
mass in the volume to the flow through the 
volume. Combustor pressure is calculated from the 
gas state Eq. (4). 

In stall operation, the fuel-air mixture in 
the combustor becomes rich. This is due to the 
decrease in airflow from the compressor. Unburnt 
fuel remains after combustion. This causes a 
decrease in combustor efficiency. The stall 
combustor model computes degraded combustor 
efficiency as a function of fuel-air ratio. 
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Turbines 


The turbine models differ significantly from 
the models of previously reported engine simu- 
lations. Previous models have used overall 
steady-state performance maps for the turbines. 
Usually, a turbine speed parameter and a turbine 
pressure ratio are used to compute turbine flow 
and enthalpy drop parameters. In the present 
model, a technique described in Ref. 11 is used. 
The turbine flow parameter, tfp, and enthalpy drop 
parameter, thp, are computed from functions of 
turbine pressure ratio only. Adjustments are made 
to account for temperature effects. The general 
flow and enthalpy drop equations for the turbines 
are as follows 

*TB = (tfp) * (11) 

Ah TB = (thp * e TB ) (12) 

Cooling bleed for each turbine is assumed to 
re-enter the flow stream at the discharge of the 
turbine. The cooling bleeds do not perform any 
turbine work. The form of the energy equation for 
the turbines requires the calculation of 
enthalpies. As previously mentioned, the 
enthalpies are calculated as functions of tempera- 
ture, neglecting the effect of fuel-air ratio. 

The form of the energy equation for the turbine 
volumes is 



Turbine pressure and temperature are 
calculated from the gas state Eq. (4). 


During stall, unburned fuel may enter the 
turbines with the possibiltiy for combustion. 
Depending on where the combustion takes place and 
on how the heat release affects the turbine work 
output, considerable model variation is possible. 
Results from various efforts to model this 
combustion process have not indicated that the 
process significantly affects simulation results. 
Because of the complexity of a turbine combustion 
model and the real-time constraint of this simu- 
lation, the turbine combustion process was not 
included in this model. 

Nozzles 

Based on available data, pressure losses in 

the bypass and the core ducts were assumed to be 
small. These pressure losses were lumped into the 
core and bypass nozzles. The nozzle flow co- 
efficients were adjusted to provide the correct 
nozzle flows for the given pressure ratios. The 
general equation form used to calculate flow for 
the bypass and core nozzles is 

*N = < cd N * A N * c % * P N )V¥ (14) 

The adjusted nozzle flow coefficient, cd'N, 
is fit as a function of nozzle pressure ratio. 

The nozzle flow parameter, cfpM , is calculated from 



(i5 > 

Rotor Dynamics 

The following general equations define the 
rotor dynamics for the high speed and low speed 
rotors. 

d/dt (N) = (30/rI) * (Q TB - Q FC ) (16) 

Q tb = (60 * J/2 * it) * ( A h TB * w tb )/N (17) 

«FC- <» * * ') * COW - ",„) F C 


* Wp c )/« (18) 

Implementation 


The engine model was programmed to run on one 
of the NASA Lewis Research Center’s two hybrid 
computer systems (fig. 3). Each system consists 
of an Electronic Associates Inc. (EAI ) PACER 600 
hybrid computer and assorted peripherals. Each 
hybrid computer consists of one PACER 100 digital 
computer, two model 681 (or 680) analog consoles, 
and a model 693 hybrid interface. A multi- 
variable function generation (MVFG) unit is 
available as a peripheral to the analog computers. 

The TF34-GE-100 turbofan engine simulation 
uses two model 681 analog consoles, one model 680 
analog console, and the MVFG unit. The consoles 
and the MVFG unit are connected to each other 
through the Central Trunking System. Figure 4 
shows the computational breakdown among the three 
consoles. Some univariate functions are imple- 
mented using digitally controlled function gener- 
ators, DCFG's that are available on the analoq 
consoles. 

To attain the real-time response of the simu- 
lation, the MVFG unit is used for bivariate 
function generation. Use of the digital computer 
for evaluation of bivariate functions required a 
time-scaling of the analog by a factor of 100. 

The digital computer introduced a delay that 
prevented stable, real-time operation of the 
simulation. Because of the use of the MVFG's, the 
digital computer is not necessary for the simu- 
lation calculations. It is used for set up and 
check out of the analog consoles and MVFG unit 
prior to simulation operation. 

The digital computer is also used to obtain 
listings of steady-state simulation data. A NASA 
developed software executive EXEC1 and an inter- 
active data collection and display program 
INFORM 1 *-* 1,5 are used in the simulation. 

Selected engine variables are read into the 
digital computer via analog-to-digital converters 
and are processed by INFORM. Guidelines for use ’ 
of the EXEC1 and INFORM software are given in Ref. 
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Concluding Remarks 

A full-range, real-time hybrid simulation of 
the TF34-6E-100 turbofan engine with post-stall 
capability is currently being developed. The 
engine model equations are based on data obtained 
from the engine manfacturer's steady-state digital 
simulation of this engine. Dynamic forms of the 
continuity, energy, and momentum equations have 
been added to provide the required transient 
capability. Where necessary, simplifications to 
the basic equations have been made to attain real- 
time response of the simulation. As of the 
writing of this paper, the out-of-stall simulation 
is running in real-time on the hybrid computer. 

The in-stall characteristics for the fan and 
compressor have been incorporated into the simu- 
lation and are being checked out. Once the simu- 
lation is fully operational, it will be validated 
by experimental data from engine stall tests that 
are planned. Eventually this simulation will be 
used to test and evaluate various stall recovery 
control schemes that are now under development. 
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Figure 1. - TF34-GE-100 turbofan engine schematic. 



(a) Normal engine airflow. 



(b) Stalled engine airflow. 
Figure 2. - Flow diagram. 
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Figure 4. - Computational breakdown of simulation. 
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